Thin front propagation in random shear flows 
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Front propagation in time dependent laminar flows is investigated in the limit of very fast reaction 
and very thin fronts, i.e. the so-called geometrical optics limit. In particular, we consider fronts 
evolving in time correlated random shear flows, modeled in terms of Ornstein-Uhlembeck processes. 
We show that the ratio between the time correlation of the flow and an intrinsic time scale of the 
reaction dynamics (the wrinkling time t^) is crucial in determining both the front propagation speed 
and the front spatial patterns. The relevance of time correlation in realistic flows is briefly discussed 
in the light of the bending phenomenon, i.e. the decrease of propagation speed observed at high 
flow intensities. 

PACS numbers: 47.70.Fw,82.40.Ck 
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Front propagation in fluid flows is a problem rele- 
vant to many areas of science and technology ranging 
from combustion technology to chemistry Q and ma- 
rine ecology In the last years several theoretical 
II H B II WU 13 and experimental [H E H H 
works studied chemically reactive substances stirred by 
laminar flows. This problem, though considerably sim- 
pler than the case of turbulent flows 1], is non trivial 
and displays a very rich and interesting phenomenology. 
We mention here the front speed locking phenomenon 
in time dependent cellular flows, which was numerically 
and theoretically found in Ref. 9] and then experimen- 
tally observed in Ref. jlSj]. Another interesting example 
is represented by the theoretical studies on time periodic 
shear flows |^ ^3 and the recent experimental work T^l 
which study aqueous reactions in periodically modulated 
Hele-Shaw flows. 

Laminar flows are interesting also because they con- 
stitute a theoretical laboratory to study some problems 
which can be encountered in more complex (turbulent) 
flows. For instance, this is the case of time correlations 
[isl that are believed to be very important in deter- 
mining the bending of turbulent premixed flame velo city 
when the intensity of turbulence is increased (see |l7| 
for a discussion about this problem). Actually for bend- 
ing several mechanisms have been proposed like reaction 
quenching dynamics of p ockets of material which 
did not react left behind [l^ and finally time correla- 
tions mm. 

The aim of this paper is to investigate the role of time 
correlations in the propagation of reactions in random 
shear flows. In particular, we shall consider the problem 
in the so-called geometrical optics, or Huygens regime 
[21I that is realized in the case of very fast reactions, 
taking place in very thin regions. As in Refs.[T5llT6|. we 
neglect possible back-effects of the transported reacting 
scalar on the velocity fleld, i.e. we treat the problem 
in the context of passive reactive transport. The latter 
assumption is justified for dilute aqueous auto-catalytic 



reactions and more in general for chemical reaction with 
low heat release. It should be also remarked that in the 
chosen framework pockets cannot be generated due to 
shear geometry and quenching of reaction cannot hap- 
pen due to the choice of working in the geometrical optics 
limit. Therefore, the case under consideration allows us 
to focus on the effects due to time correlations solely. 

Let us start to introduce our problem by shortly re- 
calling the main equations. Since we consider premixed 
reactive species, the simplest model consists in study- 
ing the dynamics of a scalar field 0{x,t) representing the 
fractional concentration of the reaction's products (6* = 1 
inert material, 9 — fresh one and < 6* < 1 coexistence 
of fresh material and products). The evolution of 9 is 
ruled by the advection-reaction-diffusion equation: 



dtU + u-V9 = DA9 



fiO) 



(1) 



where it is a given velocity field (incompressible V • u = 
through this paper). The f{9) (that is typically a non- 
linear function with one unstable ^ = and one stable 
9=1 state) models the production process occurring on 
a time-scale r^. 

Eq. 1^ may be studied for different geometries and 
boundary conditions. In this work we consider an infi- 
nite two-dimensional stripe along the x-direction with a 
reservoir of fresh material on the right, inert products on 
the left and periodic boundary conditions in the trans- 
verse direction (which has size L). In particular, we shall 
be concerned with the concentration initialized as a step, 
i.e. 9{x, y,0) = 1 for x < 0, and zero otherwise. With this 
geometry a front of inert material (stable phase) propa- 
gates from left to right with an instantaneous velocity 
which can be defined as: 



dx dt9{x, t) , 



(2) 



more precisely this is the bulk burning rate (integra- 
tion is over the entire domain V) . Most of the theoretical 
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studies aim to predict the dependence of the average front 
speed Vf = (vf) on the details of the velocity field. Very 
important are of course also the propagation speed fluc- 
tuations; one would like to predict how these are related 
to the fluid velocity fluctuations. These are in general 
very difficult issues, but very important in technological 
applications, where one has to project the reactor geom- 
etry and flow characteristics. Definite answers about the 
reaction propagation exists only in particular conditions, 
e.g. when the flow is motionless (u = 0) and under rather 
general hypothesis on the production function f{9) it is 
possible to show that the reaction asymptotically propa- 
gates with a velocity vq within the bounds (see Ref. 
for an exhaustive review): 



Dm D (f{e)\ 



(3) 



where /' indicates the derivative, and the thickness of the 
reaction zone varies as ^ oc y^/Jrv. For a wide class of re- 
action terms /, such as the autocatalytic reaction dynam- 
ics, f{9) = 9{1 — 9), and more in general for convex fimc- 
tions (/ (9) < 0) one can prove that vq — 2^J Df'{Q)/Tr 
exactly. In the presence of a velocity field it, generally 
one has that the speed Vf is larger than the bare velocity 
Vq. Specifically here we consider the limit in which the 
reaction is much faster than the other time-scales of the 
problem, formally this regime is reached when ^ 
and D ^ but with D/tj. = const such that the bare 
propagation velocity vq = jr^ is finite and well de- 
fined, while the reaction zone thickness shrinks to zero 

0) [23 , where for the sake of notation simplicity we 
posed /'(O) = 1.. It should be noted that this regime, 
also called geometrical optics limit is commonly encoun- 
tered in many applications [2H . In this limit, being sharp 

0), the front dynamics can be described in terms 
of the evolution of the surface (line in 2c?) which divides 
inert (9 = V) and fresh (9 = 0) material. The effect of 
the flow is thus to wrinkle the front increasing (in two 
dimensions) its length Cf and, as a consequence of the 
relation [1II2H 



(4) 



(where L is the length of a flat front in the absence of fluid 
motion) its propagation velocity, i.e. > vq. Quantify- 
ing such an enhancement is one of the main goals of, e.g., 
the community interested in combustion propagation [l| . 
It should be also remarked that the presence of compli- 
cated flow has also an important role in the generation 
of patterns, i.e., front spatial structures. 

From a formal point of view the evolution of 9 can be 
recast in terms of the evolution of a scalar field G(r,t), 
where the iso-line (in 2d) G'(r, i) = represents the 
front, i.e. the boundary between inert (G > 0) and fresh 
(G < 0) material. G evolves according to the so-called 
G-equation 



The analytical treatment of this equation is not trivial, 
and even in relatively simple cases (e.g., shear fiows) nu- 
merical analysis is needed. Also on the numerical side 
solving (j^l is a non trivial issue, indeed the presence 
of strong gradients usually requires the regularization of 
JSJ through the introduction of a diffusive term (see e.g. 
(23). Here, following Ref. 0], we adopt a Lagrangian 
integration scheme the basic idea of which is now briefly 
sketched. 

First of all let us introduce the type of flow we are 
interested in, we consider shear flows that can be written 



as 



u={U{t)g{y),Q), 



(6) 



being g(jj) the functional shape of the flow (here g{y) — 
sin(27ry/L)) and U{t) its intensity. The domain of inte- 
gration is chosen as a stripe [0 : A^i] x [0 : L], where 
N (typically in the range 5 — 20) is the maximum num- 
ber of cells of size L in the x-direction that are used in 
the integration (the number should be fixed according to 
the front width). The number of cells N is dynamically 
adjusted. In particular, after the propagation sets in a 
statistically stationary regime, while the front propagates 
the cells on the left that are completely inert with 6* = 1 
are eliminated by the integration domain. On the right 
side we retain only a finite number (which depends on 
the maximum allowed speed) of cells with fresh material 
9 — 0. The domain is discretized and the value of 9 in 
each point of the lattice is updated with the following 
rule. At each time step, each grid point r„^m = (xn, ym) 
is backward integrated in time according to the advection 
by the flow dr/dt = —u. Once the point r' that will ar- 
rive in Tn^m at time t is known, 9{rn^m, t) is set to 1 if in 
a circle centered in r' and having radius vodt there is at 
least one grid point with 9 — 1. This is a straightforward 
way to implement the Huygens dynamics. The algorithm 




dtG 



VG^vo\VG\ 



(5) 



FIG. 1: Typical front patterns for a stationary shear flow 
(a), time correlated shear flow with Tf — 200 (b) and with 
Tf = 2 (c). In the stationary case U — l/\/2, while in the 
time correlated case we set Urms = \/2. In all cases the bare 
velocity is vq = 0.2. For (a) and (b) we used Ny = 800 grid 
points and Ny = 3000 for (c). Here and in the following 
figures L = 2iT. 
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FIG. 2: Normalized velocity Vn = ~ vo)/U vs t/t^ 

for: vq = 1, U = 1 (circles, green online), vq = 0.2, U = 1 
(squares, red online) and vo = 0.5, U = 1 (triangles, blue 
online). The corresponding t^'s were numerically computed 
as W}/2U obtaining 2.6, 9.48 and 18.96, respectively {W} is 
estimated by counting the number of pixel in the border be- 
tween inert and fresh material) . The inset shows the unsealed 
results. The resolution used is Ny = 800. 
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FIG. 3: Measured front velocity vf{t) (thick black line) and 
the adiabatic prediction no -I- |(7(i)| (solid line, red online) 
versus t/r/ for a correlated flow with Urms = 1, iiq = 0.2 and: 
(top) Tf = 200(> t» = 31.4), (bottom) r/(2 < t„ = 31.4). 
The resolution used was Ny = 800 in the first case and Ny = 
3000 in the second one. 



works as soon as v^dt is sufficiently larger than the spatial 
discretization dx = dy = L/Ny (where Ny is the number 
of grid points in the y-direction, and N^ = N Ny). For 
a detailed description of the algorithm see the Appendix 
in Ref. 0. For a stationary shear flow, i.e. U{t) = U, 
by means of simple geometrical reasonings one can show 
that at long times the front evolves with velocity 



Vf = vo + U sup{g{y)} , 

V 



(7) 



which, with the choice of the sin flow, means Vf = vq + U. 
Similarly one can predict the asymptotic shape of the 
front, which is shown in Fig.^. The important features 
are the presence of a stationary (maximum) point in cor- 
respondence of the point where g(jj) has its maximum, 
and a cusp in its minimum. The asymptotic speed Q is 
reached only after the transient time necessary to the 
front shape to reach its maximum length (corrugation). 
Following we call as the wrinkling time, that can 
be defined as the time the front width Wf (i.e. the dis- 
tance between the leftmost point in which 9 = and the 
rightmost in which 6=1) employs to pass from the ini- 
tial zero-value (indeed at the beginning the front is flat) 
to the asymptotic one . For U ^ vq this time can be 
estimated as 



tw (X L/vq . 



(8) 



This comes from the fact that starting from the flat pro- 
file the front width W/(t) grows in time as 2Ut up to 
the moment in which the cusp (see Fig. QJi) is formed 
(see also ^). Then the growth slows down up to the 
stationary value yV|. Assuming the linear growth up to 
the end one may estimate = Wf/{2U). Further, since 
in the shear flow case (where the formation of pockets 
of inert material is not possible), for U ^ vq, the width 



is proportional to the stationary front length Cf, 
which is linked to the asymptotic velocity by Eq. Q). 
Finally since the latter given by (|7|) one ends up with 

= {L/U){1 + U/vq) which reduces to ijSJl for [/ > vq. 
In Fig. 121 we show Vf{t) as a function oit/t^, as one can 
see with this rescaling the asymptotic speed is reached at 
the same instant for systems which have different U and 
Wo, as the nice collapse of the different curves indicates 
(compare with the inset). We noticed that as soon as 
U/vq > 4 (X L/vq as predicted by Eq. 

The wrinkling time is an inner time scale of the re- 
action dynamics, which is very important when con- 
sidering time-dependent flows. In particular, here we 
study the example of random shear flows lO with g(y) = 
sin(27r?//i) (as in the stationary case) and random am- 
plitudes U{t) which are chosen according to an Ornstcin- 
Uhlembeck process. Therefore, U evolves according to 
the Langevin dynamics 



'dt 




(9) 



where ?7 is a zero-mean Gaussian white noise and Tf 
defines the flow correlation time so that {U{t)U{t')) = 
U^„^s exp(— |t— i'l/ry). Clearly one has to distinguish two 
limiting cases: i) when the flow fluctuations are slower 
than the wrinkling time: Tf 3> t^; ii) when they are 
much faster: Tf <$^tw. 

i) In this condition the front has enough time for adi- 
abatically adjust itself on the instantaneous flow veloc- 
ity. Thus by generalizing Q it is natural to expect that 
Vf{t) = Vo + \U{t)\ (as confirmed in Fig. |3Ji) so that 
Vf = Vq + {\U {t)\) . In other words if the velocity fluctu- 
ations are slower than the wrinkling time the front can 
be efficiently corrugated close to the maximal wrinkled 



4 




1 2 3 4 5 6 7 



FIG. 4: Normalized velocity Vn — {Vf — vo)/{\U\) as a func- 
tion of Tf/tw for Wo = 0.2 and Urms = 2%/2 (circles, green 
online), Urms ~ \/2 (squares, red online) and Urms ~ l/'v/S 
(triangles, blue online). The inset displays the normalized 
variance cr„ in the same cases. The resolution used goes from 
Ny — 3000 (for the lowest value of Tf) up to 800 (for the 
highest one). 

shape allowed by the flow and so by Q can reach maxi- 
mal speed. 

a) On the other hand in opposite limit Tf <C the front 
has not time to be maximally corrugated by the flow, 
and so its speed cannot reach the maximal amplification 
allowed by the fluid. In this case it is not anymore true 
that Vf{t)^vo + \U{t)\ (see Fig.|3)). 

These effects on the propagation speed have a coun- 
terpart in the patterns of the front. This is evident by 
looking at the front shape (compare Fig. ^ and c with 
a). Indeed while in the case r/ ^ tw at any instant the 
shape of front closely resembles that obtained in the sta- 
tionary case, when t/ <C one notices that the front 
length is strongly reduced and the spatial structure com- 
plicated by the presence of more than one cuspid. 

Looking at Fig. it is clear that the reactive dynam- 



ics acts as a sort of filtering of the fluid velocity so that 
not only the front speed is not enhanced at the maxi- 
mal allowed value but also its fluctuations are much de- 
creased. In Fig. 0] we show the normalized front speed 
Vn = {Vf — vo)/{\U{t)\) and the normalized variance 
an = (Tf ly/{\U{t)\^) - {\U{t)\)'^ (i.e. the standard de- 
viation of the Vf{t) normalized by that expected on the 
basis of the adiabatic process |[/(t)|), by fixing the flow 
intensity Urms and varying the correlation time. Note 
that in the limit of very long correlation times Vn ~ 1 
and cr„ « 1 . As one can see a fast drop of the front speed 
and average fluctuations with respect to its maximum 
value is observed when r/Zt^, < 1, confirming the above 
picture. 

These results along with those of Refs. 0, ITHl 
confirm the importance of time correlations in the fiow 
in determining the front speed. This may be relevant 
to more realistic flows in the light of the above cited 
bending phenomenon. Indeed in turbulent flows one has 
that at increasing the turbulence intensity fluctuations 
on smaller and smaller scales appear. These are charac- 
terized by faster and faster characteristic time scales. In 
this respect, as suggested in by the results of this work, 
one may expect that in the corrugation by these scales 
become less and less important, so that the average front 
speed may increase less than expected. 

We conclude by noticing that it would be very inter- 
esting to test the effects of time-correlations on the front 
propagation also in other kind of laminar flows. In par- 
ticular this could be performed in experimental studies 
in settings similar to those of Refs. [131 where flows of 
the form u{r, t) = U{t)v{r) can be easily generated with 
a good control of the time dependence of the amplitude. 
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